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^ ' Abstract 



We use field theoretic renormalization group methods to study the 
'^ I critical behavior of a recently proposed Langevin equation for driven 

lattice gases under infinitely fast drive. We perform an expansion 
CJ ■ around the upper critical dimension, dc = 4, and obtain the critical 

exponents to one loop order. The main features of the two loop calcu- 
lation are also outlined. The renormalized theory is shown to exhibit 
a behavior different from the standard field theory for the DLG with 



5t ! finite driving, i.e. it is not mean field like. 

Since it was first introduced by Katz et al. |]l]], the driven lattice gas model 
(DLG hereafter) has attracted considerable interest [^ ^. Being one of the 
simplest archetypes of non-equilibrium model its study may contribute to 
pave the way for an understanding of out of equilibrium systems. The DLG 
consists of a periodic regular lattice on which nearest-neighbor particle-hole 
exchanges are performed. The hopping rate is determined by the energetics 
of the Ising Hamiltonian H, the coupling to a thermal bath at temperature T, 
and an external uniform driving field E pointing along a specific lattice axis. 
In particular, the hoping rate depends on [(Aif + iE)/T], where AH is the 



energy variation which would be caused by the tried configuration change, 
E = |E|, and i = 1 (-1) for jumps along (against) E and otherwise (see 0, Q 
for a detailed description). The DLG exhibits a continuous phase transition 
from a disordered state at high T to a stripe like ordered state at sufficiently 
low T p, ^. The nature and properties of this transition have been largely 
studied in recent years. A new general Langevin equation has been proposed, 
designed to capture the physics of the DLG at the critical point |Q. This 
Langevin equation has different relevant terms for the cases < E < oo 
and E' = oo (in which particles cannot jump against the field) respectively 
[^, Q. That is, the point E = oo behaves as a sort of tricritical point in 
the parameter space where some terms are exactly zero and the relevance 
of the different operators has to be re-evaluated. For finite values of the 
driving field E the Langevin equation previously proposed by Janssen and 
Schmittmann |^ is recovered. It is the purpose of this paper to investigate 
the critical behavior of the DLG for E' = oo in order to determine explicitly 
whether the differences with the < ii^ < oo case are of a relevant nature 
and whether, therefore, the critical behavior is changed. 
The new Langevin equation reads 0, |^ 
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where Vy (V_l) is the gradient operator in the direction parallel (perpendic- 
ular) to the electric field, and the noise satisfies 

(l(x,t)) = 0, 
(V-|(x,t)-V'-|(x',t')) = -V'6{^-^')6{t-t'). (2) 

This equation is analogous to a model B in the direction perpendicular to the 
field (where the energy takes into account the interaction with the parallel 
direction through the crossed derivatives term), coupled to a simple random 
diffusion mechanism in the parallel direction. This means that the relevant 
ingredient of the infinite driving is the anisotropy it introduces, while the di- 
rectionality of the flux is irrelevant from a renormalization group point of 
view. A very similar equation has been proposed to describe the Freeder- 
icksz transition in nematic liquids, and general asymmetric two-dimensional 
pattern formation M. 
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Figure 1: Elements of perturbation theory: the response and correlation 
propagators and the four-point vertex. 0-legs are indicated by a wiggle. 



In order to renormalize this equation, following standard methods |p , let 
us introduce a Martin-Siggia-Rose response field and recast Eq. (|l|) as a 
dynamical functional [H , the associated action of which is 
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The free propagators are 
G°2(k,a;) 
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and the vertex is: —605^/12 k^. These elements can be represented diagram- 
matically as in Figure 1 (wavy legs symbolize response fields; straight lines 
stand for density fields). 

In order to renormalize the theory, one has to look for the primitive 
divergences in a perturbation expansion. If Tn^ denotes a one-particle irre- 
ducible vertex function with n external 0-legs and n external 0-legs, only Fn 
and Fi3 are found to possess primitive divergences. The Feynman diagrams 
contributing to these vertex functions (shown in Figure 2) are topologically 
identical to model B graphs [§. However, the bare correlation and response 
propagators that follow from Eq. (^) are anisotropic, in contrast to their 
counterparts in model B g]. 

To one loop in e = 4 — d, the ultraviolet divergences in Fn and F13 
lead to the renormalization of r and u, the latter being the dimensionless 





Figure 2: One loop diagrams contributing to Fn (left) and F13 (right). 

coupling constant u = AgT~'^^'^g. A^ is a numerical factor to be defined 
below. We define renormalized parameters tr and m^j by tr = Z^-t and 
Mij = ZuU. Given that the leftmost diagram in Figure 2 does not depend 
on external moments nor frequencies, therefore the derivatives of Fn with 
respect to them vanish, and no extra (field) renormalizations are required. 
The Z factors are determined by the following normalization conditions 
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A convenient choice for the normalization point NP is ki 

T = /i^, where // is an arbitrary momentum scale. To one loop, we find 

Fn(w,k) = tuj + ^kUe + t) + Di, 



(5) 
and 



eO; 



T,,{w,k) = -^kig + D,, 



(6) 



where Di{D2) corresponds to the algebraic expression of the left(right) dia- 
gram in Figure |^. A calculation in dimensional regularisation |jlO|, P] yields 
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Figure 3: Two loop contribution to F 
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The renormalization group equation obtained after requiring invariance of the 
bare irreducible vertex functions upon changes on the normalization point 
reads 



lida + I3du^ + C^. 
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where the renormalization group functions are defined in the usual way: 
I3{ur) = fidfj,UR, and C{ur) = /i(9^(lnri?). A straightforward calculation 
then leads to 
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from which one can determine the location and stability of the fixed points. 
To this order, apart from the trivial mean-field result mIj = 0, a nontrivial, 
infrared stable, fixed point -ujj = 5 emerges. This fixed point controls the 
critical behavior of the theory below four dimensions. 

Now we proceed with the calculation of the associated critical exponents. 
We first note that, as indicated above, no renormalization of the fields 0, 



has been required. Therefore, in particular, the anomalous dimension of 
vanishes up to one loop, i.e., rj = + O^e"^). Concerning the exponent 
z/^, which controls the divergence of the correlation length with temperature 
pip, we simply have z/^ = C('"r)~S and z/^ = 1/2 + e/4 + 0{e^). This is to 



be compared with z/ = 1/2, the value obtained by Janssen and Schmittmann 
in the standard field theory 0. This result demonstrates that the continuous 
version Eq. ^ of the DLG with E = oo is not mean- field like but char- 
acterizes a universality class other than the one in j^. Since there are no 
dangerous irrelevant operators in Eq. (^) standard scaling laws apply (con- 
trary to the case in 0). Therefore the exponents are related to each other 
and estimating t] and z/^ is sufficient to deduce all the other exponents. For 
instance, the order parameter exponent (3 can be written as /? = ^{d — 2-\-rj) 
|ig], and we have /? = 1/2 + 0{e^). 

The previous results concern the one loop approximation. The two loop 
calculation presents an interesting new feature, namely, that the scaling be- 
comes fully anisotropic. In fact. Figure 3 reveals that, contrary to what hap- 
pens to one loop order, the two loops correction to Fn depends on external 
frequencies and momenta; and in absence of any symmetry between paral- 
lel and perpendicular derivatives, one can easily convince himself by simple 
inspection that df^iTu ^ df^idf^iTu. In order to absorb these two different 

divergences, one is constrained to renormalize the parallel and the perpendic- 
ular momenta in a different way: if k± —>■ lk± then k\\ — >■ l^~^'^k\\ (where 7 oc £^ 
can be determined by explicitly computing the derivatives of the diagram in 
Figure 3). The scaling law has to be rewritten as /5 = h'±/2{d — 2 + 7 + 77), 
and all the exponents become non mean-field in this approximation. 

Summing up, we have performed the renormalization of the field theory 
in for the DLG under an infinitely large driving field. The renormalization 
procedure yields results essentially different from those for a finite field. In 
particular, corrections to mean field are observed explicitly in the one loop 
approximation for the exponent z/^. Anisotropic exponents and a non- mean- 
field exponent /3 appear from simple arguments based on the analysis up to 
two loop diagrams. These severe differences with respect to the finite driving 
field case call for extensive computational simulations to observe numerically 
the physical differences between both cases. 
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